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PREFACE 
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1.0  INTRODUCTION 


Increased  interest  in  the  study  of  the  fluid  mechanics  and  rate 
kinetics  of  condensation  processes  in  expansion  flow  fields  has  re¬ 
sulted  from  a  heightened  desire  both  to  prevent  their  occurrence, 
such  as  in  aerodynamic  simulation  facilities,  and  to  exploit  them  for 
such  applications  as  molecular  beam  studies  and  the  possible  use  of 
the  injection  of  charged  clusters  for  refueling  of  fusion  devices. 

Many  investigations  (Refs,  1  through  3)  have  used  mass  spectrometric- 
type  sampling  diagnostics  which  are  necessarily  located  in  the  far- 
field,  low  density  region  of  the  expansion  and,  consequently,  distant 
from  the  spatial  region  of  condensation  onset  and  growth.  Therefore, 
such  measurements  observe  the  results  of  the  integrated  rate  kinetic 
processes  along  the  sampled  streamtube.  It  is  obviously  desirable  to 
obtain  nonperturbing,  spatially  resolved  measurements  in  the  vicinity 
of  the  growth  processes  to  obviate  the  need  for  inferral  of  information 
from  downstream,  far-field  data,  and  optical  scattering  techniques 
offer  such  a  possible  diagnostic  approach. 

Although  light  scattering  was  used  as  early  as  1951  (Ref.  4)  to 
study  condensation  in  expansion  flows,  the  technique  did  not  exhibit  its 
power  until  the  early,  inadequate  light  sources  were  replaced  by  laser 
sources.  The  most  notable  application  of  laser  scattering  to  such 
condensation  studies  includes  that  of  Wegener  and  his  students,  who 
emphasi2ed  investigations  of  the  condensation  of  an  impurity  species 
in  a  non -condensing  gas,  i.  e. ,  the  isothermal  condensation  process 
(Refs.  5  through  7).  Further,  Beylich  (Refs.  8  and  9)  has  applied 
this  technique  to  the  study  of  CO2  condensation,  and  Daum  and  his 
associates  (Ref.  10)  have  studied  the  condensation  of  air  expansions. 

In  these  referenced  works  where  appropriate  data  were  available, 
the  condensate  radii  ranged  in  value  from  10  to  approximately  100  A 
with  number  densities  from  10^^  to  10^4  cm indicating  that  the 
scattering  is  adequately  described  by  the  Rayleigh  theory  rather  than 
the  more  complicated  Mie  formulation.  With  these  results  used  as  a 
basis,  laser  Rayleigh  scattering  diagnostics  have  been  used  at  AEDC 
to  study  a  variety  of  condensing  flows  (Refs.  11  through  14).  Further, 
Raman  scattering  (Refs.  13  through  16)  has  been  used  for  the  measure¬ 
ment  of  both  monomer  density  and  temperature  throughout  the  flow 
field,  including  the  region  of  the  onset  of  condensation.  The  purpose 
of  the  experiments  reported  herein  was  to  use  laser  Rayleigh  scat¬ 
tering  as  a  diagnostic  to  locate  the  axial  onset  of  condensation  in  an 
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expansion  flow  field  and  measure  the  spatial  growth  characteristics  of 
the  condensed  clusters.  With  systematic  variation  of  gas  species, 
reservoir  pressure,  orifice  or  throat  diameter,  and  nozzle  expansion 
angle,  these  measurements  provided  determination  of  the  methods  of 
scaling  condensation  phenomena  with  gas  species  and  source  reser¬ 
voir  conditions  and  geometry. 


2.0  RAYLEIGH  SCATTERING 


Since  the  basic  equations  and  rationale  for  the  application  of 
Rayleigh  scattering  to  the  study  of  condensing  gas  flow  fields  is  pre  - 
sented  in  Refs.  11  through  14,  only  a  brief  summary  is  given. 

For  an  incident  laser  beam  of  wavelength  X  and  intensity  I© 
focused  within  a  gas  sample  of  number  density  N  with  species  polar¬ 
izability  O',  the  scattered  intensity,  I,  which  is  normalized  by  Iq,  is 
given  by 


I  =  K^Na^A*  (1) 

where  the  constant  K^,  collects  all  unimportant  transmission  and  cali¬ 
bration  factors.  For  a  scatterer  of  radius  "a"  which  is  characterized 
by  bulk  properties,  it  is  known  that  a  a  a^,  which  indicates  the  sensi¬ 
tivity  of  the  process  to  scatterer  size.  If  one  assumes  the  condensing 
flow  field  to  be  composed  of  a  collection  of  gas  phase  monomers  and 
molecular  clusters,  or  i-mers,  where  i  is  the  number  of  molecules 
per  cluster,  the  single  Rayleigh  scattering  intensity  with  polarization 
parallel  to  the  incident  beam's  plane  of  polarization  is  (Ref.  11) 


r(ii)  = 


(2) 


where  Nq  is  the  reservoir  number  density.  The  scattered  intensity 
I'  (II)  includes  the  further  normalization  provided  by  the  scattered 
intensity  from  a  collection  of  monomers  of  number  density  Nq. 


For  an  uncondensed,  isentropic  expansion. 


(3) 
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Super-  and  subscript  zeros  denote  isentropic  and  reservoir  conditions, 
respectively.  The  axial  variation  of  I°(n)  is  provided  by  the  method- 
of -characteristics  solution  (MOCS)  (Ref.  17)  for  nozzle  flow  and  by  the 
Sherman -Ashkenas  theory  (Ref,  18)  for  sonic  orifice  flow. 

Deviation  of  the  measured  1*  (II)  from  I°(II)  indicates  the  existence 
of  condensation,  A  direct  measure  of  the  existence  of  clusters  within 
the  flow  is  given  by  the  scattering  function,  f,  which  is  written  as 


where  is  the  i-mer  mole  fraction.  The  approximation  results  be¬ 
cause  for  small  values  of  condensate  mass  fraction  it  is  assumed  that 
Ni  »  Nrp,  where  N-p  is  the  total  local  number  density.  Although 
the  scattering  function  (f)  is  an  ambiguous  measure  of  the  simultaneous 
increase  in  the  mean  cluster  size  and  condensate  mole  fraction,  the 
axial  variation  of  f  as  a  function  of  Pq,  D  or  D^,  di/2>  and  gas  species 
yields  empirical  condensation  scaling  laws. 

Assuming  small  clusters  to  be  characterized  by  weak  chemical 
bonding,  the  cluster  polarizability  is  assumed  to  be  additive;  i.  e. , 
orj  *  iai.  Cons’equently,  for  a  monodisperse  distribution  in  cluster 
sizes,'  i  =  J  and  aj  «  Jocj,  so  that  f  is  J^Xj.  The  condensate  mass 
fraction,  g,  is 

6  -  JX/[T  +  (J  -  DXj) 


The  depolarization  ratio,  p,  of  the  mixture  of  monomers  and 
J-mers  is  the  ratio  of  the  measured  scattered  intensity  components 
polarized  perpendicular  and  parallel  to  the  incident  plane  of  polar¬ 
ization;  that  is. 


P  =  r(i).'['(n)  (6) 

3.0  EXPERIMENTAL  APPARATUS 


Figure  1  shows  the  experimental  configuration.  A  cylindrical 
vacuum  chamber  approximately  1 . 2  m  in  diameter  and  3 . 0  m  in  length 
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which  was  equipped  with  LN2  and  20K  gaseous  He  pumping  to  provide 
a  low  pressure,  10“*^  to  lO"^  torr  background  environment  for  the 
vacuum  flow-field  expansion  studies.  The  flow  generators  were 
mounted  on  a  motor-driven,  three-dimensional  movement  mechanism 
with  an  accuracy  and  reproducibility  of  0.013  cm.  The  sonic  orifices 
were  of  1.325-,  3.2-,  and  3.05-mm  diameter  with  a  diameter-to- 
thickness  ratio  greater  than  20.  The  conical  nozzles  had  a  nominal 
throat  diameter  and  nozzle  length  of  1.0  mm  and  5.  334  mm,  respec¬ 
tively;  the  three  nozzle  half-angles  were  14.  5,  9.0,  and  5.  63  deg. 
Schematics  of  the  sonic  orifice  and  conical  nozzle  are  shown  in  Fig.  2. 
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Figure  1.  Experimental  arrangement. 


The  gas  reservoir  was  instrumented  with  standard,  calibrated 
pressure  and  temperature  gages.  Gases  were  supplied  from  high 
pressure  bottles,  and  the  stated  purity  of  the  gases  used  was  as 
follows:  N2,  99.  998  percent;  O2,  99.  5  percent;  CO,  99.  5  percent; 
Ar,  99.99  percent.  No  further  purification  was  performed;  to 
minimize  effects  of  particulate  matter,  two  25.0-nm  filters  were 
installed  in  the  inlet  line. 

As  shown  in  Fig.  1,  the  flow  is  in  the  x-direction,  the  laser 
beam  injection  is  along  the  y-axis,  and  scattered  radiation  is 
observed  in  the  z-direction.  For  the  Rayleigh  scattering  measure¬ 
ment  the  argon  ion  laser  was  operated  at  1.0  W  power  at  514.  5  nm; 
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however,  for  flow  visualization  4.0  W  total  laser  output  power  was 
used.  The  incident  laser  radiation  polarization  was  rotated  along  the 
x-direction,  expanded,  and  focused  onto  the  chamber  centerline. 

Light  scattered  from  the  focal  volume  was  collected  by  an  f/2  lens 
system,  collimated,  and  focused  onto  the  input  slit  of  a  0.  85-m  double 
grating  spectrometer.  For  the  Rayleigh  scattering  measurements, 
HN-22  Polaroid®  material  was  placed  in  the  collimated  light  path,  and 
a  polarization  scrambler  was  placed  immediately  in  front  of  the  spec¬ 
trometer  entrance  slit.  The  entrance  slit  aperture  setting,  collec¬ 
tion  optics  magnification,  and  beam  focusing  together  resulted  in 
observation  of  a  1.  5-mm-long,  50-  to  100-/i<m-diam  cylindrical  scat¬ 
tering  volume. 


Sonic  Orifice 


3/8  l<0.04D 

0  ■  1.325  to  3.2  mm 


Conical  Nozzle 


Measurements  in  Inches 
Except  as  Noted 

Roure  2.  Sonic  orifice  and  conical  nozzle  schematics. 

The  detector  was  a  thermoelectrically  cooled  EMI-9 502B  photo¬ 
multiplier,  and  the  output  was  processed  by  an  Ortec®  photon¬ 
counting  system  for  either  digital  display  or  strip  chart  recording. 

Several  laser  beam  input  apertures  were  used  to  reduce  back¬ 
ground  radiation  resulting  from  laser  plasma  light  and  forward 


II 


AEDC-TH-7a-67 


scatter  off  of  optical  components.  Laser  and  viewing  dumps  were  pro 
vided  for  further  reduction  of  background  signals,  and  all  optically 
accessible  surfaces  were  either  painted  with  a  flat  black  coating  or 
covered  with  a  black  flocking  material. 


4.0  RESULTS  AND  ANALYSIS 


4.1  RAYLEIGH  SCATTERING  RESULTS 

Figures  3  through  6  show  the  experimental  results  of  the  axial 
profiles  of  I'  (II)  for  the  N2,  CO,  O2,  and  Ar  expansions  from  sonic  . 
orifices.  Figures  7  through  9  show  the  axial  variation  of  I'  (II)  for 
N2  expansions  from  the  three  conical  nozzles  used.  Theoretical  pre¬ 
dictions  as  obtained  from  the  Sherman-Ashkenas  theory  or  the  MOCS 
are  also  shown.  It  is  observed  that  the  onset  of  condensation  is  mani¬ 
fested  by  a  dramatic  increase  of  I'  (II)  relative  to  the  iseritropic  pre¬ 
diction.  The  onset  of  condensation  moves  nearer  the  saturation  point 
as  Pq  increases,  and,  with  the  exception  of  the  region  of  discontinuity 
in  the  MOCS  nozzle  calculations,  I'  (II)  is  in  good  agreement  with  the 
calculated  values  prior  to  condensation  onset.  Furthermore,  it  can  be 
seen  that  for  the  lowest  Pq  values  the  metastable  gas  sample  supports 
a  supersaturated  state  for  approximately  30  throat  or  sonic  orifice 
diameters  before  condensing.  It  is  also  noted  that  the  massive  conden¬ 
sate  growth  for  the  nozzle  flows  is  quite  abrupt,  whereas  the  massive 
condensate  growth  region  for  the  sonic  orifice  flow  is  generally  pre¬ 
ceded  by  a  gradual  deviation  from  the  isentropic  prediction. 

Radial  profiles  of  I'  (II)  are  shown  in  Fig.  10  for  the  Pq  =  10.2 
atm,  01/2  ®  5 -deg  conical  nozzle  for  three  values  of  x.  The  data 

for  X  =  17,  35  are  interesting  in  that  two  scattering  peaks  are  symmet¬ 
rically  located  off  the  centerline.  Similar  observations  have  been 
reported  by  Beylich  (Ref.  8)  in  a  study  of  CO2  condensation  in  a  nozzle 
flow.  From  Fig.  9  it  is  seen  that  at  x  =  17.35,  condensation  has  not 
begun  on  the  axial  centerline.  However,  Fig.  10  shows  that  onset  has 
already  begun  for  r/Dt  >  0.  Figure  11  shows  photographic  observation 
of  the  Pq  "  10. 2  atm,  61/2  =  14.  5-deg  nozzle  flow  field,  and  the  nozzle, 
the  dark  isentropic  expansion  zone,  and  the  bright  onset  zone  are 
clearly  evident,  as  is  filamentary  structure  with  the  condensation 
growth  region. 
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Figure  3.  Axial  variation  of  I'dl)  for  all  N2 
reservoir  pressures  and  sonic 
orifice  diameters  investigated. 
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Figure  4.  Axial  variation  of  l'(ll)  for  all  CO 
reservoir  pressures  investigated. 
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Figure  9.  Axial  variation  of  I'dl)  for  all  Nj  Figure  10.  Radial  variation  of  I'dD  for  an  Na 

reservoir  pressures  investigated,  reservoir  pressure  of  10.2  atm  at 

8^  =  14.5  dog.  axial  positions,  0%  =  14.5 
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Figure  11.  Flow  visualization  photograph,  axial  scan, 
Ni,  P,  “  10.2  atm,  6y^  =  14.5  deg. 


Figures  12  and  13  show  the  axial  variation  of  the  scattering  func¬ 
tion,  f,  with  N2  gas  for  the  14.  5-deg  half-angle  nozzle  and  the 
D  =  1.325-mm  sonic  orifice,  respectively.  These  plots  are  typical 
of  all  the  f  versus  x  plots  obtained.  It  is  noted  that  in  these  semi- 
logarithmic  plots  the  f  values  for  a  given  Pq  and  nozzle  (or  orifice) 
form  straight  lines,  and  it  is  the  intersection  of  these  straight  lines 
with  the  x-axis  that  is  used  to  determine  the  onset  of  condensate 
growth.  These  axial  onset  locations  are  denoted  by  1^^.  The  rapid 
increase  in  f  following  onset  is  obvious,  as  are  the  orders  of  magni¬ 
tude  increase  in  f  as  the  reservoir  pressure  increases. 

Using  the  vapor  pressure  data  of  Hilsenrath,  et  al.  (Ref.  19)  and 
the  isentropic  solution  for  each  particular  flow  field  investigated,  the 
saturation  values  of  pressure  and  temperature,  Pg  and  Ts,  respec¬ 
tively,  are  obtained,  as  are  similar  values  at  condensation  onset, 
denoted  by  Pq  and  Tg.  The  isentropic  supersaturation  pressure 
ratio,  (sg)°,  is  defined  as 


(7) 
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Figure  12.  Axial  variation  of  scattering  function,  f,  for 
14.5-deg  nozzle,  N2  gas. 
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Figure  13.  Axial  variation  of  scattering  function,  f,  for  sonic 
orifice,  D  =  1.325  mm,  N2  gas. 
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and  the  isentropic  degrees  of  supercooling,  (s  g)®,  are  defined  as 


(sg)"--  (8) 

It  should  be  noted  that  the  supersaturation  ratio  defined  here  is  not 
the  normal  definition  of  the  ratio  of  the  pressure  at  onset  to  the 
equilibrium  vapor  pressure  at  onset,  The  (s^)®  values  are 

orders  of  magnitude  lower  than  the  normally  defined  supersaturation 
ratios  for  the  flow  fields  investigated  here.  These  supersaturation 
parameters  are  illustrated  in  Fig.  14,  a  diagram  of  the  expansion 
process. 


Figure  14.  Diagram  of  expansion  process. 


Figure  15  shows  the  variation  of  x^/xg*  (s0)°,  and  (sq)®  with  Pq 
for  the  three  conical  nozzles  studied.  The  approach  of  the  location  of 
condensation  onset  toward  the  saturation  point  as  Pq  increases  is 
obvious.  The  highest  value  of  supersaturation  ratio  for  a  given  Pq  is 
achieved  by  the  ~  63 -deg  nozzle  flow,  and  typically  40K  super- 

cooUng  is  possible.  Figures  16,  17,  and  18  show  the  variation  of 
xg/xg,  (sg)°,  and  (s  g)®  with  P^  for  the  sonic  orifices  studied.  The 
highest  value  of  supersaturation  ratio  for  a  given  Pq  was  achieved  by 
the  D  =  1.325-mm,  N2  expansion.  Supercooling  typically  ranged  from 
30  to  50K.  Table  1  is  a  tabulation  of  the  saturation,  condensation  on¬ 
set,  and  supersaturation  parameters  for  the  various  gases  and  sources. 
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Rffjre  15.  Variation  of  degrees  of  super¬ 
cooling,  and  supersaturation  ratio 
with  reservoir  pressure  for  the 
three  conical  nozzles  investigated. 
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Figure  16.  Variation  with  reservoir 

pressure  for  sonic  orifices. 
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Figure  17.  Supersaturation  ratio  versus 
reservoir  pressure  for  sonic 
orifices. 
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Figure  18.  Degrees  of  supercooling  versus 
reservoir  pressure  for  sonic 
orifices. 
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Table  1.  Saturation  and  Condensation  Onset  Parameters 
for  Gases  and  Sources  Studied 
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leegreesi 

Pl'Ofl w 
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0 
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3 

0 
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I 

i 
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IK/ 

■■o 

llor'l 

^8 
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•V" 

14.5 

1 

M 
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13  2 
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27.1 

5&5 

22  4 
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23  9 
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34.6 

?87 
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18.3 

56.8 

27.6 

a  40 

18  1 
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m 

5  05 

7  I 

11.3 

54.8 

31  5 
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7  4 

7  D 
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0046 

12.5 

152  2 
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285 

2  66 

7.9 

4.85 
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42.0 

0  0235 
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206 

40  1 
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83 

3  1 

48.4 
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9.0 
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22.6 
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38  7 

283 

6  05 

101 
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24  8 

a  275 

17.4 
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386 
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016 
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8L3 
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80 
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1  22 
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5.6 
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D.263 
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4.2  ANALYSIS  OF  SCALING 

Empirical  functional  relations  of  f  with  x,  x^,  and  D  or 
were  obtained  using  the  results  of  the  axial  variation  of  the  scat¬ 
tering  function  data.  These  axial  variations  were  represented  by  the 
form 


A 


X 


.bf 


(9) 


and  by  graphical  determination  it  was  found  that  ’’b"  could  be  repre¬ 
sented  by 


b 


KP. 


<10) 
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For  the  sonic  orifices  it  was  found  that  K  could  be  represented  by 


K  = 


(11) 


where  Cq,  Cj^,  and  n^  are  constants.  Further  graphical  results 
showed  that 


A 

*5  = 


(12) 


for  the  sonic  orifices,  where  C2,  03,  and  ni  are  constants.  Table  2 
is  a  tabulation  of  the  sonic  orifice  scaling  constants  determined 
where  there  was  sufficient  data.  '  For  the  conical  nozzles  it  was  de¬ 
termined  that 


K  = 


cot0,/2D;"^ 


(13) 


Table  2.  Sonic  Orifice  Scaling  Constants 


Gas 

0 

(mm) 

"0 

•=1 

N2 

1.325 

3.05 

2.94 

2.93 

0.222 

0.696 

2.02 

(0.25) 

O2 

1.325 

2.98 

0.513 

CO 

1.325 

2.87 

Ar 

3.2 

3.03 

0.49 

It  should  be  noted  that  the  nozzle  throat  diameter  variation  was 
extremely  small  and  was  a  result  only  of  imprecision  during  nozzle 
fabrication.  Therefore,  the  variation  of  K  with  Dt  was  mainly  in¬ 
ferred  from  the  sonic  orifice  results,  and  it  remains  to  be  experi¬ 
mentally  verified.  Again,  graphical  resxdts  showed  that 

wh'ere  the  variation  with  is  again  mainly  inferred.  Table  3  is  a 
tabulation  of  the  conical  nozzle  scaling  constants. 
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Table  3.  Conical  Nozzle  Scaling  Constants 


Gas 

(mini 

«U2 

{deg) 

%  "1  'l  <=3 

1.04 

145  ) 

1.00 

9,0  1 

2.98  0.525  (2)  (0.25 

LOO 

5.63  ) 

Values  of  nj  *0.5  and  cg  =  0.  25  are  used  to  show  the  scaling  of 
condensation  onset  with  Pq,  D,  or  Dt  and  61/2  in  Figs.  19  and  20  for 
the  conical  nozzles  and  sonic  orifices,  respectively.  It  is  seen  that 
xg  for  the  conical  nozzles  is  a  linear  function  of  (P^D^  cot  61/2)“^^^ 
and  for  the  sonic  orifices  xg  is  a  linear  function  of  (P§D)“l/4. 

Equations  (9)  through  (14)  may  be  used  to  write  the  scattering 
functions  for  the  sonic  orifices  and  conical  nozzles  as  follows^ 


f  = 


1  Cl 

D  1  en 


C4  b->  fn[ftA] 


(15) 


and 


f  =  (c d!"'  tan  (9 


l/2P>fn 


A  H  I  '  C  3 


=  c^b-^£n($A') 


(16) 


respectively.  For  values  of  Uq  =  3  and  ci  =  2  the  scaling  of  the  scat¬ 
tering  function  (or,  in  this  case,  of  experimental  values  of  b"l)  for  the 
various  gases  and  sources  is  shown  in  Fig.  21.  It  is  noted  that  these 
choices  of  scaling  constants  give  a  good  linear  variation  of  the  scat¬ 
tering  function  with  the  reservoir  and  gas  source  parameters  used. 

On  the  basis  of  far-field  mass  spectrometric  sampling  of  sonic 
orifice  jets  of  noble  gases,  it  has  been  shown  (Ref.  3)  that  the  12:6 
Leonard- Jones  potential  parameters  e/k,  e/o^,  and  a  are  appropriate 
scaling  parameters  for  the  source  temperature,  pressure,  and  orifice 
diameter,  respectively.  It  has  been  proposed  (Ref.  20)  that  dimer 
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concentration  for  expanding  jets  scales  as  Pq  it  also  been 

observed  that  light  scattering  in  regions  prior  to  massive  condensate 
growth  scales  as  PqD  (Ref.  12).  Furthermore,  the  condensation  onset 
locations,  Scg,  have  been  shown  to' be  a  function  of  Pq  Dgq  in  this  investi¬ 
gation,  With  this  in  mind  the  intermolecular  potential  constants  (Ref.  21) 
given  in  Table  4  have  been  used  to  plot  the  condensation  loci  of  the  flow- 
fields  studied  in  a  reduced  pres  sure -temperature  plot.  In  Fig,  22, 

]^g  been  plotted  versus  T=!'  where 


P'"  = 

T*  =  TgTf/k), 

and 

d;.,  =  IV. 

for  sonic  orifices,  while 

=  C(y)Dj 


(17) 

(18) 


(19) 


(20) 


for  conical  nozzles. 


Table  4.  Intermolecular  Potential  Constants 


Gas 

T'  =  E/k 
(K) 

e 

(ergs) 

D'  =  0 
(cm) 

P'  =  £/0^ 

(torrl 

N? 

95.05 

1. 31x10’^'^ 

3. 70x10'® 

1. 94x10^ 

0. 

117.5 

1. 62x10'^^ 

3.58x10'® 

2.65x10^ 

CO 

100.2 

1. 38x10'^^ 

3. 76x10'® 

1.94x10^ 

CO2 

205 

2. 83x10'*'* 

4. 07x10'® 

3. 15x10^ 

Ar 

119.8 

1.65x10'*^ 

3. 40x10'® 

3. 14x10^ 

ciy  =  1. 4)  is  0. 

36  (Ref.  22).  It  is  noted  in  Fig.  22  that  the  homo- 

nuclear  molecules  N2  and  O2  have  a  common  condensation  locus  for 

the  sonic  orifice 

expansions,  whereas  monatomic  Ar  and  polar  CO 
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have  widely  different  loci.  The  condensation  loci  for  the  nozzles  are 
of  identical  slope  but  are  displaced  from  each  other.  In  Fig.  23, 
has  been  plotted  versus  T’*"'  for  the  homonuclear  species  studied.  It 
is  noted  that  the  conical  nozzle  onset  locus  and  the  sonic  orifice  onset 
locus  are  very  close  to  each  other.  The  conical  nozzle  onset  locus 
may  be  represented  as  follows: 

.  T*  ’-' 

whereas  for  the  sonic  orifices. 


(21) 

(22) 


Figure  19.  Scaling  of  conical  nozzle  condensation  onset,  N2. 


■1(4 


Figure  20.  Scaling  of  sonic  orifice  condensation  onset,  N2 . 
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Figure  21.  Scattering  function  scaling. 
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Figure  22.  Condensation  onset  locus  using  reduced  values. 
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Figure  23.  Condensation  onset  locus  using  reduced  values. 


As  a  check  it  was  desired  to  derive  the  form  of  Eq.  (22)  for 
sonic  orifices  using  source  flow  and  isentropic  relations  and  the 
spatial  location  of  condensation  onset  relation  given  by  Eq.  (12). 
The  source  flow  relation  is 


^  =  c,(y) 


^  2-2y) 


(23) 


where 

Cj(y  =  1.4)  =  0..37.i  ) 

y  Ref.  22. 
c,(y  =  1.67)  =  0.281  ) 
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From  the  experimental  results. 


ft 

*5  = 


(24) 


where  a  functional  dependence  on  reservoir  temperature,  Tq,  has 
been  anticipated.  The  isentropic  relation  is 


p 


cr 


(25) 


For  y  =1.4,  Eqs.  (23),  (24),  and  (25)  show  the  condensation  on 
set  locus  to  be  of  the  form 


P  . 


(26) 


and  for  7  =  1.  67, 


P  «  <27 

which  compares  favorably  with  Eq.  (22)  for  y  =  1.4  and  with  the  re¬ 
sults  shown  in  Fig.  22  for  -y  =  1. 67. 

For  nozzle  flow  fields  the  situation  is  more  complex  since  the 
source  flow  relation  is  no  longer  valid.  A  typical  T/Tq  axial  vari¬ 
ation  is  shown  in  Fig.  24  for  the  14.  5-deg  nozzle.  For  each  of  the 
three  values  of  expansion  half-angle,  a  straight  line  was  fitted  to  the 
logarithmic  plots  of  the  MOCS -predicted  T/Tq  axial  variation.  The 
results  indicated  that 


-0.92B 


(28) 


was  a  fair  approximation  for  the  nozzles.  Using  Eqs.  (28),  (14),  and 
(25),  the  predicted  nozzle  condensation  locus  was 

p  -  (29) 


which  compares  favorably  with  Eq.  (21). 


30 


AEDC-TR-76-67 


+  e-Beam  Values  (Ref.  23) 
o  I'dl) 

A  Nj/Nq 
*  Tr/To 

□  Depolarization  Ratio,  p 
Pg  •  10. 2  atm.  N2 
l4.5-deg  Half-Angle  Conical  Nozzle 


Figure  24.  Axial  variation  of  number  density, 

temperature,  mass  fraction  (g),  and  depolarization 
ratio:  14.5-deg  conical  nozzle,  Pq  =  10.2  atm,  N2. 


4.3  CORRELATION  OF  TEMPERATURE.  RAYLEIGH  SCATTERING,  AND 
DEPOLARIZATION  RATIO  MEASUREMENTS  WITH  CONDENSATION 

Figures  24  through  30  show  the  measured  axial  variation  of  I'(II), 
N/Nq,  and  T/Tq  for  several  Pq  values  for  N2  nozzle  and  sonic  orifice 
flows.  The  number  density  and  temperatures  were  obtained  using 
laser-Raman  scattering  diagnostics  (Refs.  13  through  16)  and  also  by 
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using  electron  beam  fluorescence  diagnostics  for  the  specific  flow 
conditions  of  Fig.  24.  It  is  noted  that  the  monomer  number  density 
is  little  affected  by  the  condensation  process;  however,  the  tempera¬ 
ture  can  be  increased  by  as  much  as  50  percent  above  the  isentropic 
prediction  due  to  the  heat  release  in  the  condensation  process.  It  is 
also  observed  that  the  increase  in  temperature  correlates  very  well 
with  the  onset  of  condensate  growth. 


o  rill) 

o  Nj/No 

*  Tr/To 

Pq  ■  6. 80  atm.  N2 

14.5-deg  Half-Angle  Conical  Nonie 


Figure  25.  Axial  variation  of  l'(ll),  number  density, 
temperature,  and  mass  fraction:  14.5-cieg 
conical  nozzle,  Po  =  6.80  atm,  N2 . 
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Figure  26.  Axial  variation  of  l’(ll),  number  density.  Figure  27.  Axial  variation  of  l'(ll)  and  temperature: 

temperature,  and  mass  fraction:  14.5-deg  sonic  orifice,  D  =  1.325  mm,  Po  =  7.44  atm, 

nozzle.  Po  =  3.40  atm,  N,.  atm.  N^. 
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O  I'dll 
•  Tr'To 


D  •  1.25  mm 
Pr»-Pure  N2 
Pg*  2.79  atm.  N2 


Figure  30.  Axial  variation  of  l'(ll)  and  temperature: 

sonic  orifice,  D  =  1.325  mm,  =  2.79 
atm,  N2. 


Depolarization  measurements  of  the  Rayleigh  scattering  were 
performed  to  acquire  some  information  concerning  the  asymmetry  of 
the  scattering  clusters.  The  axial  variation  of  the  depolarization 
ratio,  p,  for  the  Pq  =  10.2  atm,  0i/2  =  14.  5-deg  conical  nozzle  con¬ 
dition  is  sho^vn  in  Fig.  24.  It  is  seen  that  p  does  indeed  decrease 
rapidly  from  its  room  temperature  monomer  value  as  one  proceeds 
axially  through  the  cluster  growth  region.  This  behavior  is  intui¬ 
tively  expected  because  as  the  linear  N2  molecules  cluster  together 
they  should  form  more  spherical  scatterers  that  contribute  a  larger 
portion  to  the  scattered  intensity  than  the  monomers.  It  is  interesting 
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to  note  that  the  depolarization  ratio  axial  variation  for  Ar  (Ref.  11) 
shows  an  initial  rapid  increase  in  p  followed  by  a  slow  decrease, 
Again  this  is  expected  because  as  the  spherical  monomeric  Ar  atoms 
cluster  to  form  dimers,  the  depolarization  ratio  should  initially  in¬ 
crease;  however,  with  further  clustering  the  depolarization  ratio 
should  decrease. 


44  COMPARISON  OF  CONDENSATION  CALCULATIONS 
WITH  MEASUREMENTS 


In  order  to  acquire  information  regarding  condensate  mass  frac¬ 
tion  and  cluster  size  along  the  centerline  of  the  expansion  flow  field, 
a  liquid-drop,  monodisperse  distribution  condensation  model  was 
used  (Refs.  23,  24,  and  25).  This  calculation  (developed  by  M, 

Kinslow)  assumes  the  condensing  flow  field  to  be  inviscid,  adiabatic, 
and  one-dimensional,  with  no  mass  transfer  across  the  streamtube 
boundary.  The  gas  and  condensate  are  assumed  to  obey  the  perfect 
gas  relation.  The  condensed  phase  is  assumed  to  be  of  the  form  of 
monodisperse  spherical  drops  or  particles  which  are  characterized 
by  bulk  properties  and  to  be  in  the  free  molecular  flow  regime  relative 
to  the  uncondensed  phase.  Condensate -gas  velocity  slip  effects  are 
ignored;  and  the  condensate  growth  rate  is  determined  by  gas  conden¬ 
sate  interaction  only.  The  mass  accommodation  coefficient  is 
assumed  to  be  unity.  The  initial  size  and  number  density  of  spontane¬ 
ous  nucleation  sites  were  adjustable  parameters  and  were  selected  to 
reproduce  as  closely  as  possible  the  experimental  results. 

The  condensation  calculation  used  the  flow-field  conditions  at  the 
saturation  point,  Xg,  as  the  starting  data.  For  the  14.  5 -deg  nozzle, 
condensation  calculations  were  made  for  reservoir  pressures  of 
10,2,  6.8,  and  3.4  atm,  and  the  results  are  shown  in  Figs,  24 
through  26,  respectively.  A  mole  fraction  of  1.  33  x  lO^S  for  dimers 
was  used  to  obtain  the  fair  agreement  between  the  calculations  and 
measurements .  These  calculations  yield  scattering  functions  within 
50  percent  of  the  measured  values.  The  measured  axial  variations 
of  number  density  and  temperature  were  obtained  using  laser-Raman 
diagnostics  (Refs.  13  through  16),  and  also  by  using  electron  beam 
fluorescence  diagnostics  for  the  specific  flow  conditions  of  Fig,  24,  as 
previously  mentioned.  It  is  seen  that  the  predicted  temperature  increase 
due  to  condensation  is  in  fair  agreement  with  the  measured  increase. 
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Also  shown  in  Figs.  24  through  26  is  the  predicted  axial  variation 
of  condensate  mass  fraction,  g.  For  Pq  =  10,2  atm,  x  =  30,  the  calcu¬ 
lations  indicate  71  molecules  per  drop  with  a  radius  of  10.  6  A.  It 
should  be  noted  that  using  these  molecules /drop  values  with  associated 
mass  fractions,  Eq.  (5)  and  the  relation  f  =  J^Xj  give  scattering  func¬ 
tions  an  order  of  magnitude  smaller  than  measured.  This  indicates 
that  the  concept  of  polarizability  additivity  is  not  applicable  for  clusters 
this  large. 


5.0  DISCUSSION  AND  SUMMARY 


With  regard  to  the  experimental  data  and  MOCS  calculations  shown 
in  Figs.  7  through  9  and  24  through  26,  it  is  noted  that  the  MOCS  shows 
a  discontinuity  in  slope  for  the  axial  profDes  of  I'{II)  or  N/N^  and  T/Tq. 
This  is  due  to  expansion  effects  from  the  nozzle  lip.  The  experimental 
data  shown  in  Figs.  24  through  26  would  indicate  that  the  discontinuity 
in  slope  is  not  realized  in  practice. 

As  noted  previously  mass  spectrometric  sampling  (Refs.  1,3, 
and  20)  has  shown  that  dimer  concentration  scales  as  P^Deq.  It  has 
been  observed  from  this  investigation  that  the  condensation  onset 
location  (more  specifically,  the  condensate  growth  onset)  scales  as 
(Po^eq^”  '  These  results  imply  that  the  termolecular  dimer  for¬ 
mation  mechanism  is  the  rate -controlling  process  for  condensation  in 
rapidly  expanding,  low  density  flow  fields.  Remembering  that  bimo- 
lecular  reactions  scale  as  Po^eq  3)  and  assuming  the  growth  of 

clusters  of  trimers  and  larger  species  to  proceed  by  bimolecular 
collisions,  one  would  expect  the  scattering  function,  f,  to  scale  as 
Po^eq*  'The  sonic  orifice  results  of  this  investigation  show  a  P^D^ 
scaling  for  the  scattering  function;  however,  the  nozzle  results  show 
a  p3(]^)/cot  ©1/2  scaling  for  the  scattering  function,  which  is  in  dis¬ 
agreement  with  the  equivalent  nozzle  concept  of  Ref.  3 .  Furthermore, 
the  condensation  onset  results  of  this  investigation  indicate  that  for  a 
given  Pq  and  Tq  a  decrease  in  61/2  with  nozzle  length  remaining  the 
same  will  result  in  decreasing  x^,  whereas  calculations  of  Sherman 
and  Griffin  (Ref.  24)  show  that  a  decrease  in  01/ 2  with  nozzle  length 
changing  to  maintain  a  constant  area  ratio  will  result  in  an  increasing 
X0.  However,  the  calculations  of  Clark  (Ref.  26)  for  constant  nozzle 
length  and  throat  diameter  show  qualitative  agreement  with  the  conden¬ 
sation  onset  results  of  this  investigation;  quantitative  agreement  is 
difficult  to  assess  because  of  the  different  definitions  of  condensation 


37 


AEDC-TR-76-67 


onset.  Future  scaling  studies  wUl  be  devoted  to  variation  of  Tq,  D^, 
and  0yj2  with  constant  area  ratio. 

With  regard  to  the  use  of  intermolecular  parameters  for  normali¬ 
zation  of  the  condensation  onset  loci  in  the  P-T  plane,  it  is  quite  inter¬ 
esting  that  the  homonuclear  N2  and  O2  show  a  common  locus'while  the 
monatomic  and  polar  diatomic  species  show  different  loci.  This  re¬ 
sult  is  in  agreement  with  the  results  of  Yealland,  et  al.  (Ref.  2), 
although  Hagena  and  Obert  (Ref.  3)  suggest  that  the  method  would  be 
restricted  to  monatomic  species.  Future,  investigation  of  this  normal¬ 
ization  procedure  will  involve  NO,  HCl,  NH3,  and  H2O. 

Further  development  of  condensation  calculations  for  flow  fields 
is  needed  in  order  to  completely  predict  the  scattering  efficacy  for 
expansions  of  a  variety  of  gases  throughout  the  flow  field.  Although 
the  present  model  makes  a  fair  prediction  of  the  axial  variation  of 
light  scattering,  it  requires  adjustment  of  initial  nuclei  size  and  mole 
fraction  to  fit  experimental  data.  Furthermore,  the  assumption  of 
unity  accommodation  coefficients  and  constant  nuclei  mole  fraction  is 
doubtful.  The  classical  capillarity  theory  of  condensation  can  be  used 
to  calculate  nucleation  rates  and  hence  the  size  and  number. of  "critical 
nuclei"  at  any  point  in  the  flow  field.  Shown  in  Fig.  31  are  nucleation 
rates  calculated  for  three  Pq  values  for  the  14.  5 -deg  nozzle  using  the 
classical  nucleation  theory.  It  is  noted  that  nucleation  has  passed 
through  a  maximum  and  rapidly  declined  prior  to  observed  condensate 
growth,  indicating  the  sequential  nature  of  condensation  processes  in 
rapid  expansions,  that  is,  saturation  followed  by  nucleation  followed 
by  condensate  growth.  However,  the  capillarity  theory  is  not  beyond 
criticism  in  spite  of  some  of  its  successes,  because  the  theory  is 
thermodynamic  in  nature  and  not  kinetic.  Ftirthermore,  variation 
of  surface  tension  values  with  radius  at  small  nucleus  sizes,  extremely 
large  super  saturation  ratios  such  as  are  foimd  in  rapidly  expanding 
flows,  and  statistical  mechanical  correction  factors  can  all  cause  great 
variation  in  calculated  nucleation  rates  (Ref.  -  23).  An  alternate  approach 
has  been  described  by  Dorfeld  and  Hudson  (Refs.  27  and  28).  Their 
method  is  strictly  kinetic  in  nature,  and  it  provides  a  method  for  calcu¬ 
lating  dimer  concentration  using  termolecular  rate  equations  and  also 
accounts  for  the  effective  increase  in  the  termolecular  collision  rate 
at  low  temperature  due  to  the  presence  of  loosely  bound  orbiting  pairs 
(Ref.  27).  Condensate  growth  is  calculated  using  bimolecular  collision 
rates  of  monomers  with  rate -limiting  species  such  as  dimers  (Ref.  28). 
It  is  this  approach  that  will  be  explored  for  future  condensation  calcu¬ 
lation  development. 
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Figure  31.  Classical  nucleation  rates  for  three  N2  expansions, 
14.5-deg  conical  nozzle. 
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NOMENCLATURE 

Badius  of  scatter er 
Badius  of  i-mer 

A  function  of  reservoir  pressure  and  orifice  diameter 
and  expansion  half  angle 

Scaling  constants 

y  -dependent  parameters  whose  numerical  value  is 
given  in  the  text 

Sonic  orifice  diameter 

Conical  nozzle  throat  diameter 

Equivalent  diameter  defined  as  D  for  sonic  orifices 
and  c(y)  Dt  cot  0i/2  for  conical  nozzles 

D  normalized  by  the  intermolecular  potential  param¬ 
eter  CT 

Bayleigh  scattering  function  defined  by  Eq,  (4) 
Condensate  mass  fraction 

Relative  Rayleigh  scattering  intensity  defined  by 
Eq.  (1) 

Relative  Rayleigh  scattering  intensity,  normalized  to 
the  relative  Rayleigh  scattering  intensity  of  a  gas 
sample  of  number  density  Nq,  polarized  parallel  and 
perpendicular  to  the  plane  of  polarization  of  the  inci¬ 
dent  lasher  beam,  respectively 

Number  of  molecules  per  cluster 

Constant  in  Eq.  (1) 

Method-of-characteristics  solution 

Number  density  of  gas  species 

i-mer  number  density 

Reservoir  number  density 

Isentropic  monomer  number  density 

Total  local  number  density 

Scaling  constants 
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Po 

Ps 

Pg 

Pye 

r 

(S0)° 

To 

Tb 

T0 

T* 

t.p. 

Xi 


A 

X0 

a 

Q-i 

Y 

e/k,e/cr^ 

h/2 

X 

p 


Reservoir  pressure,  atm 
Saturation  pressure 

Pressure  at  condensation  onset 

Equilibrium  vapor  at  condensation  onset 

Pg  normalized  by  the  intermolecular  potential 
parameter  e  /  <7^ 

Radial  distance  from  flow-field  centerline 

Isentropic  super  saturation  ratio 

Isentropic  degrees  of  supercooling 

Reservoir  temperature,  K 

Saturation  temperature 

Temperature  at  condensation  onset 

Tg  normalized  by  the  intermolecular  potential  param¬ 
eter  e/k 

Thermodynamic  triple  point 
i-mer  mole  fraction 

Axial  position  in  flowfield,  x  =  x/D  for  sonic  orifices 
and  X  =  x-f;/  Dt  for  conical  nozzles 

Axial  location  of  saturation 

Axial  location  of  condensation  onset 

Polarizability 

i-mer  polarizability 

Specific  heat  ratio 

Intermolecular  potential  parameters  given  in  Table  4 
Expansion  half -angle  of  conical  nozzle 
Wavelength  of  scattered  radiation 
Depolarization  ratio  of  Rayleigh  scattering 
Intermolecular  potential  parameter  given  in  Table  4 


■44 


